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O . ABSTRACT 
(N 

Aims. Bright HNC 1-0 emission, rivalling that of HCN 1-0, has been found towards several Seyfert galaxies. This is unexpected 
O , since traditionally HNC is a tracer of cold (10 K) gas, and the molecular gas of luminous galaxies like Seyferts is thought to have 

bulk kinetic temperatures surpassing 50 K. There are four possible explanations for the bright HNC: (a) Large masses of hidden cold 
gas; (b) chemistry dominated by ion-neutral reactions; (c) chemistry dominated by X-ray radiation; and (d) HNC enhanced through 
\l ' mid-IR pumping. In this work we aim to distinguish the cause of the bright HNC and to model the physical conditions of the HNC 

and HCN emitting gas. 

Methods. We have used SEST, JCMT and IRAM 30m telescopes to observe HNC 3-2 and HCN 3-2 line emission in a selection of 
5 HNC-luminous Seyfert galaxies. We estimate and discuss the excitation conditions of HCN and HNC in NGC 1068, NGC 3079, 
NGC 2623 and NGC 7469, based on the observed 3-2/1-0 line intensity ratios. We also observed CN 1-0 and 2-1 emission and 
discuss its role in photon and X-ray dominated regions. 
Q . Results. HNC 3-2 was detected in 3 galaxies (NGC 3079, NGC 1068 and NGC 2623). Not detected in NGC 7469. HCN 3-2 was 

• detected in NGC 3079, NGC 1068 and NGC 1365, it was not detected in NGC 2623. The HCN 3-2/1-0 ratio is lower than 0.3 only 

^ I in NGC 3079, whereas the HNC 3-2/1-0 ratio is larger than 0.3 only in NGC 2623. The HCN/HNC I-O and 3-2 line ratios are larger 

, than unity in all the galaxies. The HCN/HNC 3-2 line ratio is lower than unity only in NGC 2623, which makes it comparable to 

galaxies like Arp 220, Mrk 231 and NGC 4418. 

Conclusions. We conclude that in three of the galaxies the HNC emissions emerge from gas of densities n < 10^ cm"^, where the 
chemistry is dominated by ion-neutral reactions. The line shapes observed in NGC 1365 and NGC 3079 show that these galaxies have 
^ , no circumnuclear disk. In NGC 1068 the emission of HNC emerges from lower (< 10^^ cm"^^) density gas than HCN (> 10^ cm"^^). 

. Instead, we conclude that the emissions of HNC and HCN emerge from the same gas in NGC 3079. The observed HCN/HNC and 

' CN/HCN line ratios favor a PDR scenario, rather than an XDR one, which is consistent with previous indications of a starburst 

' component in the central regions of these galaxies. However, the A'(HNC)/A'(HCN) column density ratios obtained for NGC 3079 

I can be found only in XDR environments. 
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>< ; extragalactic tracer of molecular gas ^densities ^(Hz) ^r th^^atio mngis from'o.03 to 0.4 i"n w7rm cot^^^^ 

a ■ 101^"' -g ™ can be as high as 4.4 in cold cores (T, < 15 K). 

■ ■ ■ ' Kohno I2005I I. The HCN to CO intensity ratio varies signifi- „, , , j- ,^ , , ■ , r 

, r ; ,,.rv- 1 1 • J • 1 1 The CN (cyanogen radical) molecule is another tracer or 

cantly, rrom 1/3 to 1/40 in starburst galaxies, and It has not been , , ^ . r c\ i j 

J J , ,- •• J j ij 1 dense gas, with a lower (by a factor or 5) critical density than 



1. Introduction based only on this scheme would predict then an HNC/HCN 

^ , , ., , , . , , ratios; 1 . However, the observed HNC/HCN abundance ratios 

The hydrogene cyanide, HCN molecule, is commonly used as an significantly between different kinds of molecular clouds 



determined whether this variation depends on the dense molec 
ular gas content or on the abundance and/or excitation condi 



HCN. CN is also chemically linked to HCN and HNC by pho- 



, j,.^. , ^ ■ J- 1 iT^xT todissociations (e.g. Hirota & Yamamoto 1999). Surveys of the 

tions. In addition, recent results seems to indicate that HCN may , rv . r^xr j h^t^ u i. j • j \ . 

, i j .cij 11 .-^1-0 transition of CN and HNC have been done in order to trace 

not be an unbiased tracer of the dense molecular gas content in , , j,^, • i • i • / a i. . i 

T TT-.^ 1 TTT Tr.^ ■ 1 ^/^o^ T ■ 1 a coW, dcnsc phasc of thc gas in luminous ga^xics (Aalto ct al. 

LlRGs and ULlRGs (Gracia-Carpio et al. 20061. It IS essential, J . hxt^i r^^ ■ • n j 

, „ 111 1 iT^xT ■ J 120021) . It was found that the HNC 1-0 luminosities often rivalled 

therefore, to use other molecular tracers than HCN, in order to ~7 m^xri « t-u i. . . j- ..u j r 

J J 1 1 • 1 J- • • 1 J those of HCN 1-0. These results seem to contradict the idea of 

understand the physical conditions in the dense gas. ^ c-rv • i ■ i • / 

„ . . , TT^NT warm (Ti, > 50 K) gas in the centers of luminous galaxies (e.g. 

A molecule of particular interest, for comparison with HCN, ,,r.,, , i ,„r.-i L ^-ifn i iir.r.ok i i ... 

• • iiM/-^ T-u J . f - . . n iixT/-^ . Wild et al. 119921 Wall et al. 11993ft whose IR luminosities were 

IS Its isomer HNC. The detection of interstellar HNC supports ^ ,' r. „ ~ — -r ^. , , .^xt 

, , r J ■ ■ 111- • J 1 1 suggested to originate from star formation rather than AGN ac- 

the theory of dominant ton-molecule chemistry in dark molecu- (jyjj (Solomon et al I 1992I ) 

lar clouds. Both species are thought to be created by the same ,. ' . , ... „ , , 

dissociative recombination of HCNH\ This ion can produce , ^° observations in the vicinity of the hot core 

HCN and HNC, with approximately equal abundances. Models °f ^L' f P!""i^"'.lf ^^^^'^^^ steady state and 
shock models, the HNC/HCN ratio decreases as the tempera- 

Send offprint requests to: J.R Perez-Beaupuits ture and density increase (e.g. Schilke et al. 119921 Talbi et al. 

* Research supported by the Swedish Institute as a scholarship 1996; Tachikawa et al. 2003). If a bright HNC 1-0 transition 

granted to J.R Perez-Beaupuits to pursue his M.Sc. degree at Chalmers line is nevertheless detected under these conditions, it could 

Tekniska Hogskola, Sweden be due to the following possible explanations: (a) the presence 
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Table 1. Sample of galaxies ^. 
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Galaxy 


Seyfert 


RA 


DEC 


V.,v.! 




Distance *" 




Sis (HCN) 






[hh mm ss] 


r ' "] 


[km 


'] 


[Mpc] 


["'-] 


[" 


NGC 3079 


2 


10 01 57.805 


+55 40 47.20 


1116+ 


:1 


15.0+1.1 


15x7.5 


5x5 


NGC 1068 


2 


02 42 40.711 


-00 00 47.81 


1137+ 


:3 


15.3+1.1 


30x 30 


10 X 10 


NGC 2623 


2 


08 38 24.090 


+25 45 16.80 


5549+ 


:1 


74.6+5.4 


8x8 


2.6 X 2.6 


NGC 1365 


1.8 


03 33 36.371 


-36 08 25.45 


1636+ 


:1 


22.0+1.6 


50x50 


16.5 X 16.5 


NGC 7469 


1.2 


23 03 15.623 


+08 52 26.39 


4892+ 


:2 


65.8+4.8 


8x8 


4x6 



a) The Seyfert classification, positions (in equatorial J2000 coordinates) and heliocentric radial velocities were taken from NED. 

b) The distances were calculated using the Hubble constant (Hq ~ 74.37 kms"' Mpc"') estimated by Ngeow and Kanbur (2006). 

c) The source sizes of the CO 1~0 transition line were estimated from the maps presented in (Koda et al. i200 2) for NGC 3079, (Schinnerer 
et al. 20001 for NGC 1068, (Bryant et al. 179991 for NGC 2623, (Sandgvist fT^ for NGC 1365, and (Papadopoulos & Allen 2000) for 
NGC 7469. The source sizes for the 7 = 2-1 transition line were assumed equal to that of the 7=1-0 line. For NGC 7469, the source 
size of the CO 2-1 emission estimated from the corresponding map presented by Davies et al. ( I2004t agrees well with the source size 
estimated for the CO 1-0 line. 

d) Source sizes of HCN 1-0 were estimated from the corresponding maps published in (Kohno et al. 12000b for NGC 3079, (Kohno et al. 12001 1 
and Heifer & Blitz 1995 ) for NGC 1068, (Davies et al. lOOT) for NGC 7469. The source sizes of NGC 2623 and NGC 1365 were estimated 
using proportions found in NGC 1068 (read text in §3.5). Because of their chemical link, the source sizes of the CN and HNC molecules 
were assumed the same as that of HCN. Due to the lack of high resolution maps, the source sizes corresponding to the emission of the 
higher transition lines were assumed equal to that of the 7=1-0 line. 



of large masses of hidden cold gas and dust at high densities 
(n > 10^ cm"^); (b) chemistry dominated by ion-molecule re- 
actions with HCNH^ at low density (« » 10** cm"-') in regions 
where the temperature dependence of the HNC abundance be- 
comes weaker; (c) enhancement by mid-IR pumping, also in 
low density regions where the lines would not be coUisionally 
excited; and (d) the influence of UV-rays in Photon Dominated 
Regions (PDRs) and/or X-rays in X-ray Dominated Regions 
(XDRs) at densities n > 10'^ cnT^ and at total column densi- 
ties A^H > 3 X 10^' cm"2 (Meijerink & Spaans 2005]). 

In the case of CN, observations of its emission towards the 
Orion A molecular complex (Rodriguez-Franco et al. 1998) sug- 
gest that this molecule is also enhanced in PDRs, but particularly 
in XDRs, where a CN/HCN abundance ratio larger than unity is 
expected (e.g., Lepp & Dalgarno 1 19961 and Meijerink, Spaans, 
Israel |2007] ). 

We have observed low and high transition lines of the HCN, 
HNC and CN molecules in a group of Seyfert galaxies, which 
are supposed to host both sources of power, AGN and starburst 
activity, in their central region. Our interest is to assess the ex- 
citation conditions of HCN and HNC, distinguish between the 
above possible causes of the bright HNC, and to explore the re- 
lation between the CN emission, XDRs and dense PDRs in these 
sources. 

In §2 we describe the observations. The results (spectral 
lines, line intensities and line ratios) are presented in §3. The 
interpretation of line shapes and gas distribution in the most rel- 
evant cases, as well as the possible explanations for the bright 
HNC and the modelling of the excitation conditions of HCN and 
HNC are discussed in §4. The conclusions and final remarks of 
this work are presented in §5. 

2. Observations 

We have used the James Clerk Maxwell Telescope (JCMT) in 
2005 to observe the HNC 7=3-2 (271 GHz) and the HCN 
7=3-2 (267 GHz) lines towards a sample of Seyfert galaxies. 
Observations of CN and HNC 7=1-0 (90 GHz) were made in 
2002 using the Swedish ESO Southern Telescope (SEST). HCN 
7=1-0 (88 GHz) data from literature were also used. The system 
temperatures ranged between 350 K and 430 K. In the case of 



Table 2. Beamsizes & efficiencies. 



Transition 


V 


HPBW 




Telescope 




[GHz] 


["] 






HCN 1-0 


88.632 


44 28 


0.65 0.77 


OSO IRAM 


HNC 1-0 


90.663 


55 27 


0.74 0.77 


SEST IRAM 


CN 1-0 


113.491 


45 22 


0.71 0.74 


SEST IRAM 


CO 1-0 


115.271 


44 21 


0.70 0.74 


SEST IRAM 


CN 2-1 


226.875 


23 11 


0.51 0.53 


SEST IRAM 


CO 2-1 


230.538 


20 11 


0.66 0.52 


JCMT IRAM 


HCN 3-2 


265.886 


18 9 


0.60 0.45 


JCMT IRAM 


HNC 3-2 


271.981 


18 9 


0.60 0.44 


JCMT IRAM 



a) The columns are divided in two sub-columns {left: OSO, SEST, 
JCMT, and right IRAM) indicating the telescope used to obtain the 
corresponding parameter. The IRAM 30m telescope was used to 
observe NGC 3079. The other sources were observed with the OSO 
20m, SEST and JCMT telescopes. 



CN 7=1-0 (226 GHz) the weather conditions were not so good, 
making the system temperature range between 490 K and 760 K. 
Pointing was checked regularly on SiO masers and the rms was 
found to be about 3". 

For the SEST observations we alternated between a 500 
MHz and 1 GHz backend, depending on weather Simultaneous 
observations with the 1 and 3 mm receiver were taking place. In 
addition, the data of the Seyfert galaxy NGC 3079 obtained dur- 
ing the IRAM survey in 2006 is included in this work. We used 
the software package XS (written by P. Bergman) to reduce the 
data and fit the gaussians. Beamsizes and efficiencies are shown 
in Table |2] 

The sample consists of five Seyfert galaxies of which 
two are Seyfert 1 type and three are considered mainly 
Seyfert 2. Table [T] lists the coordinates of the center po- 
sitions observed in these galaxies, their sub-classification 
as Seyfert galaxies and heliocentric radial velocities, ac- 
cording to the NASA/IPAC Extragalactic Database (NED) 
(http://nedwww.ipac.caltech.edu/). The distances were calcu- 
lated assuming a Hubble constant of Hq = 74.37 ± 2.27 
km s"' Mpc"' (Ngeow and Kanbur. 12006] ) . 
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3. Results 

3.1. NGC 1068 

The molecular line emissions observed in NGC 1068 are shown 
in Figure [1] A first-order polynomial was used in most cases to 
correct the baselines, with the exception of the HCN 3-2 and 
CN 2-1 spectra, for which a second-order polynomial was re- 
quired. The velocity resolution was set to 15 km s ' for CO and 
CN, whereas a 25 km s"'resolution was used for HNC and HCN. 
These velocity resolutions represent less than 10% of the line 
widths. The spectra are centered with respect to the heliocentric 
systemic velocity v^ys = 1 137 km s"' (from NED). 

The J-l-Q and 7=2-1 fines of CO show a triple structure 
where the peaks can be attributed mainly to the bars and the 
spiral arms, as described by Heifer & Blitz 1 19951 

We detect two of the main spingroups of CN 1-0: the 1-0 {J 
= 3/2 - 1/2, F - 5/2 - 3/2) line at the center of the spectrum, and 
the 1-0 (7 = 1/2 - 1/2, F = 3/2 - 3/2) fine shifted 856 kms"' to 
the right. This spectrum shows three components as well, how- 
ever they are hard to distinguish due to the blending. In fact if we 
freely fit three gaussian components, the uncertainties of the cen- 
ter velocity, amplitude and line width are about 100% or larger. 
Instead, if we set the line widtfi of tfie central component to 154 
kms ' , wfiicfi corresponds to tfie line widtfi found for tfie CN 
7=2-1 line as described below, we get reasonable values. Only 
tfie first gaussian component sfiows fiigfi uncertainties in tfie am- 
plitude and line widtfi. Since tfie beam sizes of tfie CN and CO 
molecules are similar at tfie frequencies of tfie 7=1-0 lines, we 
tfiink tfiat botfi beams pick up emission coming from tfie spi- 
ral arms and bars. Altfiougfi, tfie nuclear region seems to be tfie 
predominant component in tfie case of CN. 



Table 3. NGC 1068 line parameters. 



Transition 


Gaussian 


V 


n 


AV 




Component 


[km s"'] 


[mK] 


[km s"'] 


CO 1-0 


1 


1007± 1 


225.2±14.9 


58± 3 




2 


1100± 3 


265.1± 3.8 


156+11 




3 


1234± 2 


215.3±10.8 


96± 4 


CO 2-1 


1 


1007± 2 


206.8±32.5 


63± 9 




2 


1121± 4 


529.3± 8.8 


179±16 




3 


1245± 3 


222.7±37.8 


81±10 


CN 1-0 


1 


1001±23 


12.8±13.2 


94±42 




2 


1091±68 


19.4± 4.7 


154" 




3 


1213±32 


15.6± 9.5 


136±54 




4 


1946±26 


8.3'' 


154" 


CN 2-1 


1 


1093± 9 


10.4± 0.9 


154±66 




2 


1289±14 


6.2± 0.9 


154" 


HCN 3-2 


1 


1103± 7 


48.7± 24 


275+17 


HNC 1-0 


1 


1073±13 


9.7± 1.8 


114+30 




2 


1268±20 


6.7± 1.5 


170+62 


HNC 3-2 


1 


1071± 9 


12.9± 1.8 


134+23 




2 


1250±24 


3.7± 2.5 


70+58 



" The line widths were set to the value found in the main component 
of the CN 2-1 line. 

'' In order to avoid the effect of the artefact in the backend, the inten- 
sity was locked to the value expected for this spingroup. 



Due to tfie second spingroup is corrupted by noise in tfie 
backend, we set tfie amplitude of tfie gaussian to 8.3 mK, wfiicfi 
corresponds to tfie expected factor of about 0.43 times tfie ampli- 
tude of tfie main spingroup, according to tfie National Institute of 
Standards and Tecfinology (NIST). Tfie resulting central veloc- 
ity of tfie second component was 1946 kms '(~ 113.169 GHz), 
i.e., about 22 MHz sfiifted from tfie expected frequency for tfiis 
spingroup. 

In tfie CN 2-1 line, tfie two spingroups at 226.8746 GHz (2- 
1,7 = 5/2- 3/2, F = 1/2- 5/2) and 226.6596 GHz (2-1, 7 = 
3/2 - 1/2, F = 5/2 - 3/2) are severely blended since tfie sfiift is 
~300 km s ' . In order to identify tfie two spingroups we first fit 
two gaussian components to get tfie line widtfi of tfie main spin- 
group, wfiicfi is 154 + 66. We tfien set tfiis value to botfi gaussian 
components and fit again tfie otfier parameters. Tfie second spin- 
group is, in tfie optically tfiin limit, a factor 0.54 weaker tfian tfie 
main spingroup. We get a factor ~ 0.6 between tfie intensities 
obtained from tfie gaussian fit. On tfie otfier fiand, tfie resulting 
center velocity of tfie second component is 1289 km s wfiicfi 
corresponds to a sfiift of about 68 MHz witfi respect to tfie ex- 
pected frequency of tfie second spingroup. Besides tfie noise in 
tfie data, tfiis sfiift may also be produced by tfie influence of tfie 
unresolved spingroups, (2-1, 7 = 3/2 - 1/2, F = 1/2 - 1/2) and 
(2-1, J = 3/2- 1/2, F =3/2- 1/2), located in between tfie two 
main spingroups. Togetfier, tfiese inter-spingroups would pro- 
duce an intensity comparable to tfiat of tfie second spingroup. 

Tfie line sfiape and intensity of tfie CN 2-1 spectrum dif- 
fers from tfie one obtained by Usero et al.'2004'(tfiereafter U04). 
Besides tfie lack of baseline coverage observed in tfie U04 spec- 
trum, tfiere is a substantial discrepancy between tfie estimated 
source sizes. In U04 tfie emission of most of tfie molecules was 
assumed to emerge from a 6"x4" region, wfiicfi corresponds to 
tfieir estimate of tfie size of tfie circumnuclear disk (CND), based 



Table 4. NGC 1365 line parameters. 



Transition 


Gaussian 


V 


n 


AV 




Component 


[kms"'] 


[mK] 


[kms"'] 


CO 1-0 


1 


1532+ 2 


322.7+5.2 


136+4 




2 


1709+ 3 


260.3+5.0 


142+5 


CN 1-0 


1 


1511+13 


16.8+3.5 


130+25 




2 


1685+33 


8.9+14 


212+92 




3 


2356+18 


4.7+1.2 


130" 


CN 2-1 


1 


1526+11 


10.4+1.5 


130" 




2 


1727+35 


8.4+1.9 


212'' 




3 


1846+35 


4.11+3.5 


130" 


HCN 3-2 


1 


1534+ 5 


51.4+34 


136+13 




2 


1706+25 


11.4+3.3 


142'^ 


HNC 1-0 


1 


1539+11 


16.0+1.7 


142+23 




2 


1750+34 


5.1+14 


195+98 



" The two spingroups are supposed to have the same line width, so 
these were set according to the value found for the main spingroup 
of the CN 1-0 line. 

'' The line width of the second component of the double peak struc- 
ture was set to the corresponding value found in the CN 1-0 line. 
The second component of HCN 3-2 seems to be affected by noise, 
so its line width was locked at AV = 142 kms"', i.e., the corre- 
sponding line width of the CO 1-0 line. 
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Table 5. NGC 3079 Hne parameters. 



Transition Gaussian V T* AV 

Component [kms '] [mK] [kms"'] 



CO 1-0 


1 


932± 4 


129.4+11.8 


143+ 9 




2 


1146± 3 


374.3+ 3.3 


269+13 




3 


1312± 6 


145.1+13.1 


74+ 9 




4 


1389± 5 


186.3+ 9.1 


72+ 5 


CO 2-1 


1 


964+15 


233.3+18.9 


175+18 




2 


1062±18 


224.3+84.9 


98+18 




3 


1160±20 


321.7+35.3 


126+20 




4 


1311±14 


428.3+ 9.4 


190+ 7 


CN 1-0 


1 


967± 9 


5.6+ 0.6 


162+23 




2 


1287±15 


4.5+ 0.6 


215+33 


CN 2-1 


1 


997+14 


8.1+ 1.5 


137+34 




2 


1251+25 


5.3+ 1.3 


203+69 


HCN 1-0 


1 


1005+30 


9.7+ 1.4 


249+73 




2 


1274+21 


10.7+ 1.7 


188+49 


HCN 3-2 


1 


1033+42 


6.4+ 2.1 


176+93 




2 


1287+32 


10.4+ 1.8 


219+70 


HNC 1-0 


1 


983+17 


6.5+ 1.7 


129+41 




2 


1236+21 


7.8+ 1.5 


157+38 


HNC 3-2 


1 


996+48 


7.4+ 3.2 


178+124 




2 


1197+43 


7.5+ 3.8 


120+90 



on the CO 2-1 high resolution map presented by Schinnerer et 
al. (l2000l l. 

If we correct the peak antenna temperature (T^ » 9.7 mK) of 
our CN 2-1 spectrum, in order to obtain the main beam bright- 
ness temperature r^b in the same way done in U04 (r„,fc - 
r* X Tj^Ji^ X D.,„h X Q^'), we get a peak T^b of about 0.42 K, 
which is about 0.18 K (or a factor of ~1.8) larger than the peak 
temperature obtained by U04. 

Due to the chemical link between HCN and CN, here we 
rather estimate the source size based on the high resolution map 
of HCN 1-0, pubUshed in Heifer & Blitz ( fl^5] l. We estimate 
that the emissions of HCN, CN (and very likely HNC as well) 
emerge from a more extended region of about 10"xlO". This 
source size can also be inferred from the HCN 1-0 map pre- 
sented by Kohno et al. (1200 11 1. 

By correcting the peak antenna temperature (as described in 
§3.5) we get a peak T^b of about 0.12 K, which is a factor 2 
smaller than what was found by U04. 

The HCN 3-2 line has its center velocity at V =s 1103 
kms ' which coincides with the main HCN emission shown in 
the position-velocity (p-v) map by Tacconi et al. ( 119941 ) (this is 
discussed in §4). The line width and line shape of this spectrum 
are consistent with HCN 1-0 spectra published in previous work 
(Nguyen-Q-Rieu et al.TW CuiTan et al.'2000^. 

The HNC spectra present a double peak profile in which the 
main peak has its center velocity at V ^ 1070 kms"' in both 
the J-l-0 and J -3-2 transitions. However, the line shape and 
intensity of the HNC 1-0 spectrum differs from the one obtained 
by Hiittemeister et al. (1 19951 1. The HNC 3-2 spectrum seems 
to be affected by noise, since the amplitude and line width of 
the secondary component have uncertainties larger than 50%. 
The different line shapes observed between the HNC and HCN 



Table 6. Line parameters of NGC 2623 and NGC 7469. The 
third column indicate the gaussian components (G.C.) used to fit 
the spectral lines. 



Source 


Transition 


G.C. 


V 




n 


AV 








[km s 




[mK] 


[kms"'] 


N2623 


HNC 3-2 


1 


5465+ 


:34 


3.4+ 0.7 


337+88 


N7469 


CO 2-1 


1 


4773: 


to 


0.15+ 


126+0 






2 


4927: 


to 


0.16+ 


123+0 



spectra would indicate that they emerge from regions of different 
kinematics. The line parameters are summarized in Table [3] 

3.2. NGC 1365 

Figure |2] shows the molecular line emissions observed in 
NGC 1365. Without considering the second spingroups of CN, 
all the molecular spectra of NGC 1365 present a double peak 
structure, irrespective of the beam size or line transition. 

This double peak structure can also be seen in the CO 2-1 
and CO 3-2 spectra showed in Sandqvist et al. (119951) and in 
Sandqvist ( 119991 ). respectively. 

The two strongest spingroups of the CN 1-0 line are also de- 
tected in this galaxy, with a double peak structure at the nuclear 
region. The strongest peak was used as reference for the main 
spingroup. We first fit two gaussian components to the main 
group and then we set the line width of the second spingroup, 
observed at around 2400 km s"', to the value found for the main 
spingroup. The proportion found between the amplitude of the 
spingroups is about 0.3 instead of the expected value of 0.43, 
according to NIST. The data around the second spingroup are 
strongly affected by noise, so its central velocity is shifted by 
about 50 km s"' from the expected value of about 800 km s"' . 

In the CN 2-1 line, the second spingroup overlaps with the 
double structure itself. The line widths of the two spingroups are 
set according to the value found for the CN 1-0 line. Given their 
combined intrinsic line strengths, the integrated line intensity ra- 
tio of the two spingroups is expected to be about 1.85 (assum- 
ing optically thin emission). If we also set the line width of the 
second peak of the double structure at 212 km s"' , as found for 
CN 1-0, we get a proportion of about 0.4 between the two spin- 
groups. Instead, if we set the amplitude of the second spingroup 
to 5.7 mK (the expected proportion) and we fit the line width of 
the double peak structure, we get a value of 1 8 1 km s"' , which is 
about 30 kms 'less than the width found in the CN 1-0 line. In 
any case, the center velocity of the second spingroup is shifted 
by about 30 km s"' more than the expected velocity. 

The HNC 1-0 and HCN 3-2 spectra have a structure similar 
to that of the CN 1-0 spectrum. Note that in CO 1-0 the double- 
peak is more pronounced than in the CN, HCN or HNC spectra, 
where the low velocity peak dominates. The line parameters are 
summarized in Table H] 

3.3. NGC 3079 

The molecular line emissions observed in NGC 3079 are shown 
in Figure[3] The spectra are centered with respect to the heliocen- 
tric systemic velocity Vjys = 1116 kms"'. The velocity resolu- 
tion was set to 20 km s" ' for HNC and HCN, and to 1 km s" ' for 
CO and CN. Four structures are observed in the CO lines, which 
are not present in the other spectra. The CO lines are picking up 
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Fig. 1. Molecular line emissions in NGC 1068. The velocity resolution was set to 25 km s^'for HNC and HCN, and to 15 km s^'for CO and CN. 
The spectra are centered with respect to the heliocentric systemic velocity Vsys = 1137 kms"'. Emission from the spiral arms are detected in the 
CO 1-0 line. The CO 2-1 line is dominated by the emission coming from the CND. The two main spingroups of CN are detected in both J=l-0 
and J=2-l transitions. In the CN 1-0 the second spingroup is corrupted by noise in the spectrum. The two spingroups are blended in the CN 2-1 
line. The different line shapes (profiles) of the HCN and HNC spectra seem to indicate that their emissions emerge from different regions. 



extended lower density gas, compared to the CN lines (which 
have similar beam size). Hence the difference in line shape. 

The CN 1-0 line present a double peak structure. Only the 
main spingroup (J = 3/2 - 1/2, F = 5/2 - 3/2) of the CN 1-0 hne is 



observed since the second spingroup falls beyond the bandwidth 
of the backend. In the CN 2-1 line, the two main spingroups, (J 
= 5/2 - 3/2, F = 7/2 - 5/2) and (J = 3/2 - 1/2, F = 5/2 - 3/2), 
are detected. But they are severely blended and merged with the 
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Fig. 2. Molecular line emissions in NGC 1365. The spectra are centered with respect to the heliocentric systemic velocity Vsy, - 1636 kms 
The velocity resolution was set to 25 kms" 'for HNC and HCN, and to 15 kms" 'for CO and CN. A double peak line shape is observed in all 
the spectra. This structure is related with the double peak emission coming from the center of the galaxy. In the CN spectra both spingroups are 
detected. The second spingroup of the CN 2-1 line overlaps the double structure itself. 



double peak structure. In this case, only two gaussians were fit- 
ted because a third (middle) component presented large uncer- 
tainties in both the amplitude and line width, due to the noise in 
the spectrum. The two components are separated by only 254 
kms"', and the proportion between their amplitudes is about 
0.65. The difference between these and the expected values is 
attributed to the noise and the blending of the emission of the 
second spingroup. 

The HCN spectrum of the 7=1-0 and J -3-2 transitions 
agree with the HCN 1-0 spectrum obtained by Nguyen-Q-Rieu 
et al. (fT992l i. 

The HNC 1-0 spectrum is different from the tentative detec- 
tion presented by Hiittemeister et al. (1995), which is broader 
and more intense than the spectrum showed here. Besides the 
similar line shapes, our HCN and HNC spectra extends from 900 
kms"' to 1300 kms"'. Instead, the HNC spectrum showed in 
Hiittemeister et al. (1995) peaks at around 1300 kms"'. Hence, 
we believe that our HNC detections are correct. In this case, be- 



cause of the similar line shapes and widths (although the HNC 
3-2 line is strongly affected by noise, and its peak intensities 
have uncertainties between 40% and 50%.), the HCN and HNC 
emissions are likely emerging from the same gas. The line pa- 
rameters are summarized in Table |5] 

3.4. NGC 2623 and NGC 7469 

The HNC 3-2 line emission observed in NGC 2623 is shown 
in the left top panel of Figure |4] The velocity resolution was 
set to 25 kms"'. The spectrum is centered with respect to the 
heliocentric systemic velocity Vsys = 5535 kms"' (NED). After 
a total integration time of 2 hours we do not detect HCN 3-2 
emission in this galaxy. 

The right top panel of Figure |4] shows the CO 2-1 spectrum 
observed in NGC 7469. The velocity resolution was set to 15 
kms"'. The spectrum is centered with respect to the systemic 
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Fig. 3. Molecular line emissions in NGC 3079. The velocity resolution was set to 20 km s^'for HNC and HCN, and to 10 km s^'for CO and CN. 
The spectra are centered with respect to the heliocentric systemic velocity Vjys = 1116 kms"'. Emission from four structures are observed in the 
CO lines. Only the main spingroup is observed in the CN 1-0 line. The second spingroup is on the right edge of the spectrum, out of the bandwidth. 
Instead, both spingroups are observed in the CN 2~1 line, although the second spingroup overlaps the double structure of the CN emission. The 
both transitions of HCN and HNC have similar line shapes, which indicates that their emissions emerge from the same region. 



velocity Vsy, - 4892 kms" . The CO spectrum shows a double do not detect any emission. The bottom panels show the spectra 
peak structure, like the one observed in NGC 1365. We observed of the not detected lines. The line parameters of the detected 
the HNC 3-2 line for about 1 hour of integration time, but we transition lines are summarized in Table|6] 
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Fig. 4. Molecular line emissions in NGC 2623 (left) and NGC 7469 {right). The velocity resolution was set to 25 kms^'for HNC and to 15 
km s^'for CO. The spectra are centered with respect to the heliocentric systemic velocities Vsy, = 5535 and Vsy, = 4892 km s^'for NGC 2623 and 
NGC 7469, respectively. We do not detect HCN 3-2 emission in NGC 2623 nor HNC 3-2 in NGC 7469. Their observed intensities were less than 

2(7. 



3.5. Line intensities and ratios 

The velocity integrated intensities are showed in Table |7] In or- 
der to compute the line intensity ratios were corrected for the 
different beam sizes obtained with different frequencies and tele- 
scopes, according to the correction factors defined for compact 
sources in Rohlfs and Wilson (2003). In the case of extended 
sources (Q.s > i^mt) the ratios were corrected for the main beam 
filling factor, which was approximated as /,„/, - Q.s /{ils + ^mb)- 
The beams, as well as the source structures, were assumed to be 
gaussians. 

The source sizes reported in Table [1] were estimated from 
high resolution maps available in the literature. For NGC 3079 a 
source size of 15"x7.5" was estimated for the CO 1-0 emission, 
considering only intensities above 15% of the peak integrated 
intensity of the contour map presented by Koda et al. (120021 ) . 
In NGC 1068 most of the CO 1-0 emission emerges from the 
two spiral arms, with the largest extension of ~40" (e.g. Heifer 
& Blitz .1995. Schinnerer et al. ,2000) . Considering intensities 
above 20% of the peak emission of the high resolution map by 
Schinnerer et al. ("2000^ the source size of the CO 1 -0 emission in 
NGC 1068 was estimated as 30" x 30". For NGC 2623 a source 
size of 8" X 8" was estimated from the CO 1-0 map presented in 
Bryant et al. d 19991 1, which agrees well with the estimate made 
by CasoU et al. (fT98^ . 

In the case of NGC 1365 the CO emission is concentrated 
in the nuclear and bar regions (Sandqvist et al. I1995I I. Hence 
the source size of the CO emission was estimated as 50" x 50", 
corresponding to intensities above 20% of the peak emission in 
the CO 3-2 map by Sandqvist (j 1999( 1. 



From Papadopoulos & Allen ( I2000I I and Davies et al. (|2004| l 
the source size of the CO emission in NGC 7469 was estimated 
as 8"x8", which correspond to intensities above 40% of the peak 
emission in the high resolution map by Davies et al. (2004). The 
criteria of selection of the source size of CO varies depending on 
the gradient of the emission observed in the different sources. 



The source size of the high density tracers CN, HCN and 
HNC were estimated through HCN maps available in the liter- 
ature. In all these maps, we observed that the HCN 1-0 emis- 
sions emerge mainly from the nuclear region of the galaxies. 
Since there are no published maps of CN nor HNC emissions 
for the galaxies studied here, the corresponding source sizes 
were considered equal to that of the HCN emission, due to their 
chemical link. Since there are no HCN maps for NGC 1365 and 
NGC 2623, the factor VOsCHCN)/ VQs(CO) ~ 0.33 found in 
NGC 1068 was used to estimate the source size of the HCN 1-0 
emission, based on their corresponding ^^(CO). The estimated 
Q.S (HCN) are shown in Table [1] 



The line ratios were computed assuming an error of 10% 
in the reported beam efficiencies rj^t, 5% of error in the main 
beam 6„h (Table |2]i, and a 10% error in the estimated source 
sizes 0s (Table [TJ. The obtained ratios are shown in Table [8] 
From these ratios we can conclude that CO, as well as the high 
density tracers - CN, HCN and HNC - are subthermally excited. 
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Table 7. Velocity-integrated intensities ^. 



Galaxy 


/(CO) 1-0 


/(CO) 2-1 


/(CN) 1-0 


/(CN) 2-1 


/(HCN) 1-0'' 


/(HCN) 3-2 


/(HNC) 1-0 


/(HNC) 3-2 


NGC 3079 


205.8±20.6 


375.4±37.6 


2.7±0.4 


4.5±0.7 


5.7±0.8 


8.0+1.4 


2.9+0.5 


5.0+1.9 


NGC 1068 


113.7±11.4 


203.6±20.5 


9.4±1.0 


6.1±0.7 


10.0±1.7 


22.7+2.5 


3.2+0.5 


3.5+0.6 


NGC 2623 


30.2±3.2'^ 




<1.3'^ 






<0.5 


0.9+0.3 


1.8+0.4 


NGC 1365 


122.3±12.3 


181.0±18.1 ^ 


5.9±0.7 


7.5±0.9 


6.0+0.1 


15.0+2.0 


4.7+0.6 




NGC 7469 


25.6+3.8 ■= 


61.5±6.4 


1.8+0.3"^ 




1.5+0.3 




1.1+0.2'^ 


<0.8 



a) The values refer to the main-beam brightness temperature, /n,b, in [K km s"']. The errors and upper limits correspond to Icr (defined by the 
r.m.s. in the spectra) and added in quadrature to the 10% of error considered for the main beam efficiencies, 77,^1,, reported in Tableland 
by the Onsala Space Observatory (OSO). 

b) HCN 1-0 integrated intensities reported by Curran et al. (f2000ll corrected for a main-beam efficiency of 0.65, according to the on-line 
values reported by OSO. Our own data are reported for NGC 3079. 

c) Integrated intensities reported by Aalto et al. ( 20021 . The values were rounded to one decimal figure and corrected by the main beam 
efficiencies reported by OSO, 0.64, 0.45 and 0.43, of the 7=1-0 transition of HNC, CN and CO, respectively. 

d) CO 2-1 integrated intensity obtained by Sandqvist et al. (1995}. A 10% of error was assumed for the reported intensity. 

e) CO 1-0 integrated intensity derived from Curran et al. ( I2000t considering a main beam efficiency of 0.43 according to the values reported 
by OSO. 



4. Discussion 

4.1. The distribution of dense gas 

We compare the spectral shape of the high density tracers with 
HCN and CO position-velocity maps available in the literature 
in order to address the location of the dense gas. 

4.1.1. NGC 1068 

In NGC 1068 the two outer peaks of the CO 1-0 spectrum 
(Figure [T]) coincide with the maximum double peak emission 
seen in the CO position-velocity (p-v) map obtained by Heifer 
& Blitz (II995I 1. so they can be attributed to the emission emerg- 
ing from the spiral arms. The center is attributed to the emission 
emerging from the CND, as can be infeiTed from the CO 2-1 
spectrum (Figure [T]!. This center peak coincides with the max- 
imum HCN emission at 1100 km s 'seen in the HCN p-v map 
obtained by Tacconi et al., ( |1994| l. This p-v map is shown in the 
left panel of Figure|5]for comparison, along with the correspond- 
ing scaled spectra observed in NGC 1068. 

The main spingroup of the CN 1-0 line (Figure ^ shows 
a shape similar to that of the CO 1-0 line. So contributions of 
emission coming from the CND, as well as from the spiral arms. 



Table 8. Line intensity ratios. 



Galaxy 


CO ^ 


CN — 


HCNf^ 


HNCH 


NGC 3079 


0.77+0.20 


0.48+0.15 


0.18+0.06 


0.25+0.12 


NGC 1068 


0.81+0.21" 


0.19+0.05 


0.47+0.14 


0.15+0.05 


NGC 2623 


b 






0.36+0.15 


NGC 1365 


0.84+0.12'^ 


0.44+0.12 


0.42+0.12 




NGC 7469 


0.97+0.28 






< 0.20 



a) For NGC 1068 we corrected the JCMT CO 7=2-1 observations 
to the beam size of the SEST CO 7=1-0 line, applying the factor 
20/44, the ratio between the respective beams. With this, the beam 
dilution effect is cancelled in both transition lines. 

b) See Casoli et al. ( |T988] l. 

c) For NGC 1365 we also corrected the CO 7=2-1 observations to the 
beam size of the CO 7=1-0 line, applying the factor 25/44 between 
the respective SEST beams. The ratio given above is slightly higher 
than the one reported by Sandqvist et al. il995b . 



can also be inferred. The main spingroup of the CN 2-1 spec- 
trum also coincides with the maximum HCN emission around 
llOOkms 

The HCN spectra seem to contain two components, although 
we fit only one gaussian to the spectrum. The main component 
fits the region in the p-v where the strongest HCN emission is 
coming from. The secondary component might be attributed to 
the secondary peak emission of HCN, observed around 1035 
km s ' . Another smaller peak is seen in the p-v map around 1255 
km s but this component is not detected in the HCN spectra. 

In the HNC spectra, instead, two clear components are ob- 
served. The main component has a center velocity of 1073 
kms ' which lies in between the two main peaks of the HCN 
emission observed in the corresponding p-v map. With respect 
to the CO p-v map (Heifer & Blitz 1995), the main HNC com- 
ponent comes from a region where the CO emission is faint. 
Instead, the secondary component of HNC seems to emerge 
from a region around 1250 kms which corresponds to the 
secondary peak emission in the CO p-v map, and roughly to the 
third peak of the HCN emission in the respective map. Although, 
the latter is uncertain due to the noise in the spectra. 

4.1.2. NGC 1365 

On top of the middle panel of Figure |5] we show the CO 3-2 
p-v map of the central region of NGC 1365, from Sandqvist 
(II999I 1. The spectra below correspond to the re-scaled spectra 
shown in Figure |2] The double peak structure of the spectra co- 
incides fairly well with the double peak emission observed in 
the p-v map. The left components of the spectra seem to emerge 
from a region around 1530 km s"' . 

The right components of most of the spectra could emerge 
from the region around 1710 kms ', with the exception of the 
CN 1-0 line. The right component of CN 1-0 seems to emerge 
from around 1750 km s ' which, in turn, coincides with the max- 
imum emission level of the right peak in the CO 3-2 p-v map. 

The two peaks in the CO 3-2 p-v map are of about the same 
intensity, which is reflected in the '^CO 1-0 spectrum, where 
both components have intensities ~ 0.3 K. The '^CO 3- 
2 and '^CO 1-0 data obtained by Sandqvist (1999) also have 
this feature. However, the CN, HCN and HNC spectra, exhibit 
a gradient between the intensities of their high and low velocity 
components. This can also be seen in the HCN 1-0 and HCO^ 
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Fig. 5. {Top-left panel) Position-velocity (p-v) map of HCN 1-0 emission in NGC 1068, adapted from Tacconi et al. ( |1994t . The spectra below 
are the respective emission lines of Figure [T] re-rescaled to fit the velocity scale of the HCN p-v map. (Top-middle panel) Position- velocity map 
of CO 3-2 emission in NGC 1365, adapted from Sandqvist ('1999). The spectra below are the re-scaled versions of the ones shown in Figure[2] 
(Top-right panel) Position-velocity map of CO 1-0 emission in NGC 3079, adapted from Koda et al. ( I2002) . The spectra below are the re-scaled 
versions of the ones shown in Figure[3] 
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1-0 spectra obtained by Sandqvist (I1999I I. This intensity gradi- 
ent could be due to either a larger abundance of the species, or 
to a higher excitation in the corresponding region of the lower 
velocities. 

Since the observed double peak structure does not change 
with the beam size (one single feature), we think that the nucleus 
of NGC 1365 lacks a circumnuclear disk, i.e., it is consistent 
with a Seyfert 1 nucleus. 

4.1.3. NGC 3079 

The four-peak structure observed in NGC 3079 can be identi- 
fied in the CO 1-0 p-v map of Koda et al. (2002). The top-right 
panel of Figure |5] shows an adaptation of the Fig. 10 in Koda 
et al. (2002) with 0" offset along the minor axis of NGC 3079. 
Since this galaxy is almost edge-on, most of the CO emission 
is probably coming from the spiral arms, besides the nuclear re- 
gion. Instead, the high density tracers are expected to be mostly 
nuclear, as in the case of NGC 1068. This can explain the double- 
peak structure of the line shape of the CN, HCN and HNC 
molecules, in contrast to the CO lines. 

Note that in this case the CN spectra are slightly dominated 
by the lower-velocity peak, as in the case of NGC 1068 and 
NGC 1365, whereas the HCN and HNC spectra of NGC 3079 
are dominated by the higher- velocity peaks. The lower and 
higher velocity peaks are not perfectly aligned. The maximum 
separation between the peaks is ~ 80 km s"' for the lower- 
velocity peaks and ~ 60 kms^'for the higher-velocity peaks. 
The lower- velocity peaks are centered around 1000 km s 'and 
the higher- velocity peaks around 1270 kms which means the 
high density tracers tend to avoid the peaks of the CO 2-1 emis- 
sion. 

Since the line shapes of CN, HCN and HNC are similar, and 
they do not change with the beam size, their emissions likely 
emerge from the nuclear region in this galaxy, in contrast to the 
case of NGC 1068, where the line shape of the HCN emission 
changes with the beam size. 

4.2. The HCN/HNC line ratios 

NGC 1068: The HCN/HNC 3-2 & 1-0 line intensity ratios 
increase towards the CND. The ratio varies from ~2.0 for the 
lower transition lines, to ~6.5 for the higher transition. The 
large 7=1-0 beam picks up emission from both the CND and 
the starburst ring, whereas the 7=3-2 beam picks up emission 
coming mainly from the CND. This could be interpreted either 
as that the abundance ratio differs between the starburst ring 
and the CND, or that the abundance is actually the same but the 
physical conditions are different in these two regions. On the 
other hand, there could also be optical depth effects since the 
difference in the ratios is consistent with a larger optical depth in 
the 7=1-0 transition line of HNC. In section 4.3.1 we estimate 
the excitation conditions of HCN and HNC, from which we 
can derive the corresponding optical depths. As described in 
Figure |6l the optical depth of the 7=1-0 line of HNC is larger 
(starting from t=0.01) than that of the 7=3-2 line (which starts 
from T=0.003) in almost all the possible excitation conditions. 
In the case of HCN, the situation is the opposite. The lowest 
optical depth of the HCN 7=1-0 line is 0.03, whereas in the 
7=3-2 Une is 0.32. In any case, X[HCN]/X[HNC] is at least 
6.5 in the CND. According to Meijerink et al. (2007), this ratio 
can be found in gas of density rtu ~ lO^cm"^ and with PDR 
conditions, at a distance from the source of ~ 10'^ cm, if the 



Habing flux Go is about 10 , or at a slightly larger distance 
of 2 X 10'^ cm, if Go = 10\ On the other hand, in an XDR 
environment, this ratio would be found at a distance of 7 x 10'^ 
cm, if the radiation flux Ffuv is 16 erg cm"^ s"', or at a much 
larger distance of about 2 X 10'^ cm, if Fff/v = 160 erg cm"^ s"'. 

NGC 3079: The HCN/HNC ratio decreases for the higher 
transitions. The HNC 3-2 emission rivals that of HCN, making 
the HCN/HNC 3-2 line ratio only 1.6. As described above, a 
similar analysis of the expected distribution of these ratios and 
molecules, in a PDR and XDR environments, can be done based 
on Meijerink et al. (I2007I I. 

NGC 2623: In this galaxy we do not detect HCN 3-2. 
We estimate an upper limit of 0.26 for the HCN/HNC ratio 
of the 7=3-2 line. We observe that the ratio does not only 
decreases for the higher transitions, but also this is the only 
galaxy in our sample where /(HCN) < /(HNC) in the 7=3-2 
line, which makes it comparable to galaxies like Arp 220, 
Mrk 2 31 and NGC 4418, according to the recent work by Aalto 
et al. (I2007I I. They propose that the overluminous HNC can 
be explained by a pumping effect due to mid-IR background 
radiation with brightness temperatures Tb > 50 K and densities 
below critical, or due to the ISM chemistry being affected by 
X-rays. According to SchiUce et al. ( il992| ), shocks are also 
possible sources of explanation. 

The HCN/HNC 1-0 ratio is about 2.0 in NGC 1068 and 
NGC 3079, and ~ 1 .4 in the others. This indicates a brighter HNC 
emission in NGC 2623, NGC 1365 and NGC 7469. Luminous 
HNC in galaxies may have the following plausible explanations: 

a) Large masses of hidden cold gas and dust. 

If the HCN and HNC emission is emerging from gas of den- 
sities n(H2) > 10^ cm"^ then the HNC chemistry would be 
dominated by reactions like HNC -H O -> CO + NH which 
would destroy HNC at higher temperatures. Thus, at high 
gas densities, a bright HNC line would imply a considerable 
amount of cold (T^ < 24 K) dense gas. 

b) Chemistry dominated by ion-neutral reactions. 

If, however, the bulk of the HCN and HNC emission is 
emerging from gas of densities ^lO'* cm"-* then the relative 
HNC abundance may be substantial, despite the high tem- 
perature. The reason for this is that, at lower densities, reac- 
tions with HCNH+ (HCN and HNC reacts with H3 to form 
HCNH^) become more important. The ion abundance is 
higher and once HCN and HNC become protonated, HCNH^ 
will recombine to produce either HCN or HNC with 50% 
probability. This scenario is interesting since the electron and 
ion abundance is likely higher in PDRs Photon Dominated 
Regions (PDRs) (Tielens & HoUenbach, 1985). Therefore, 
in a PDR chemistry, the connection between HNC and ki- 
netic temperature may be weak since we expect the HCNH^ 
reactions to be important. 

c) Chemistry dominated by hard X-rays. 

The X-ray irradiation of molecular gas leads to a so called X- 
ray dominated region (XDR) (e.g. Maloney et al. 1996) simi- 
lar to PDRs associated with bright UV sources. The more en- 
ergetic (1-100 keV) X-ray photons penetrate large columns 
(10^^-10^"* cm"^) of gas and lead to a different ion-molecule 
chemistry. Models of XDRs by Meijerink & Spaans (2005) 
and Meijerink et al. ( 120071) indicate that the HNC/HCN col- 
umn density ratio is elevated (and larger than unity) com- 
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Table 9. Intensity ratios between the high density tracers. 



Galaxy 


HCN J O 
HNC ^ " 


HCN , ^ 
HNC -'"^ 


™- 1.0 

HNC ' " 


CN 2-1 
HNC 3-2 


OL 1-0 

HCN ^ " 


CN 2-1 
HCN 3-2 


NGC 3079 


2.15±0.67 


1.60+0.75 


0.64+0.20 


1.25+0.58 


0.30+0.09 


0.78+0.25 


NGC 1068 


2.01+0.65 


6.48+1.95 


1.98+0.58 


2.58+0.78 


0.98+0.29 


0.40+0.29 


NGC 2623 


1.4" 


<0.26 


<0.77" 




<0.5'' 




NGC 1365 


1.35±0.37 




0.87+0.25 




0.64+0.17 


0.67+0.19 


NGC 7469'' 


1.50+0.57 




0.85+0.28 




0.57+0.21 





a) Ratios reported by Aalto et al. (2002 ). 

b) New ratios computed using the HCN 1-0 source size estimated from (Davies et al. l2004t . 



pared to PDRs and quiescent cloud regions for gas densities 
around 10^ cm"-^. 
d) HNC enhanced through mid-IR pumping. 

Both HCN and HNC may be pumped by an intense mid-IR 
radiation field boosting the emission also from low density 
regions where the lines would not be collisionally excited. 
For HNC the coupling to the field is even stronger than for 
HCN, thus increasing the probability for IR pumping in ex- 
treme galaxies, such as Mrk 231. Ultraluminous galaxies, 
such as Mrk 23 1 and Arp 220, have central mid-IR sources 
with optically thick radiation temperatures well in excess 
of those necessary to pump the HNC molecule (Soifer et 
al. 1999). Even if the HNC abundance is lower than HCN, 
the HNC emission may have a higher filling factor due to the 
IR pumping (i.e. IR pumped emission from gas clouds oth- 
erwise at too low density to excite the HNC molecule) (e.g., 
Aalto et al. 2002). 



The study carried out so far allows us to distinguish between 
the above scenarios in some of the sources presented here. The 
number density n(H2) required for alternative (b) is too low to 
efficiently excite the HNC (or HCN) 1-0 and 3-2 lines, and we 
should therefore expect subthermal HCN and HNC excitation in 
this case. An HNC 3-2/1-0 line ratio of 0.3 or less is an indi- 
cation that the gas densities are below 10^ cm"-', depending on 
temperature and column density, as it is discussed in section 4.3. 
A HNC 3-2/1-0 ratio lower than 0.3 is observed in NGC 3079, 
NGC 1068 and NGC 7469 (Table®. 

The ratio ~0.4 found in NGC 2623 is in the limit between 
case (a) and (b). On the other hand, case (c) and (d) cannot eas- 
ily be ruled out. More information is required in order to distin- 
guish between the proposed scenarios, as suggested by Aalto et 
al. 120071 

The HCN 3-2/1-0 line ratio is below 0.3 only in NGC 3079, 
whereas it is larger than 0.4 in NGC 1068 and NGC 1365. 
This result is interesting since it implies that in NGC 3079 both 
HCN and HNC emission emerge from the same gas, whereas 
in NGC 1068 the HNC emission has to emerge from a lower 
(< 10^ cm"^) density gas than HCN. It would be interesting to 
see if this result holds for NGC 1365. We expect to obtain the 
HNC 3-2 data for this galaxy in a future project. 

The pumping scenario (c) should lead to a HNC 3-2/1-0 line 
ratio close to unity. This is not observed either in the HNC or in 
the HCN data we have. However, it is not possible to rule out 
this scenario since low excitation may be the result of mid-IR 
pumping of low density gas. Detailed modelling is needed in 
this case. 



4.3. Excitation conditions of HCN and HNC 

We used the radiative transfer code RADEXQ (Van der Tak et 
al. ' 20071 1 to explore a wide range of possible excitation condi- 
tions that can lead to the observed line ratios. This code is sen- 
sible to the column density of a molecule per line width, and 
uses a constant temperature and density of the collision partner, 
which in this case is H2. Another limitation of RADEX is that 
it cannot handle large optical depths (t < 100). Our analysis 
is not depth dependent and we assume a homogeneous sphere 
for the escape probability approach. Hence, our models aim to 
reproduce a sort of average cloud that represent the physical con- 
ditions of the emitting gas, which is a well fitted starting model 
for single dish observations, where all the emissions detected are 
convolved with the telescope beams. 

The grid consists of densities between 10"* and 10^ cm"^, 
temperatures between 4 and 200 K, and column densities per 
line width between 10'" and 10"^ cm"^ km"' s. Excitation maps 
were generated in order to obtain the total column density per 
line width (N/ Av cm"^ km" ' s) as function of the kinetic tem- 
perature (Tii K) and the number density of molecular hydrogen 
(n(H2) cm"^). The countour lines of the maps describe the di- 
chotomy between temperature and density. This means that, for 
a given column density, the observed 3-2/1-0 line ratio can be 
obtained with high temperatures and low densities, or low tem- 
peratures and high densities. 

The hyperfine structure of the HCN 7=1-0 transition line 
is not included in this study. The extrapolated HCN data of the 
LAMDA0 database were used instead (Schoier et al. 2005 ). Note 
that the hyperfine components may be overlapping and may in- 
teract radiatively in AGN-like environments. This process re- 
quires further analysis and modelling that is not included in 
RADEX. We are not able to generate an excitation map of CN 
due to the lack of collision data for this molecule. Although, an 
study to extrapolate the collision data for CN from other know 
molecules (e.g. CS) is ongoing. 

4.3.1. NGC 1068 

Figure |6] shows the excitation condition maps for the HCN and 
HNC molecules, modeled from the average ratios observed in 
NGC 1068. These maps show that HCN (top) requires higher 
density gas and higher column density than HNC (middle) in 
order to obtain the observed 3-2/1-0 line ratios. 

Considering a kinetic temperature Ti^ = 80 K and a density 
n(H2) ~ 10^ cm""*, as estimated by Tacconi et al. (il994 ), we 
found that A^(HCN)/Ai; ~ 1 .5 x lO''* cm"^ km"' s, which agrees 
with the respective value found by Tacconi et al. 



' http://www.sron.rug.nl/~vdtak/radex/index.shtml 
^ http://www.strw.leidenuniv.nl/~moldata/ 
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Fig. 6. Excitation conditions modelled for the 3-2/1-0 line ratios 
of HCN (top) and HNC (middle) observed in NGC 1068. The 
conditions required for the HCN and HNC molecules overlap 
in a narrow region. The relative column densities in the overlap 
zone of the excitation conditions of these molecules is shown in 
the bottom plot. The optical depth in the whole region explored 
for HCN ranges between 0.03 and 10 in the 7=1-0 line, and be- 
tween 0.32 and 30 in the J-3-2 line. In the case of HNC the 
optical depth ranges between 0.01 and 30 in the J-l-0 line, and 
between 0.003 and 30 in the the J-3-2 line. The optically thin 
limit of both molecules and lines is depicted by the right edge of 
the excitation conditions, whereas the optically thick limit cor- 
respond to the left edge of the figures above. 





1 ■■! ,i 



Fig. 7. Excitation conditions modelled for the 3-2/1-0 line ra- 
tios of HCN (top) and HNC (middle) observed in NGC 3079. 
The excitation conditions required for these molecules overlap 
in most of the range explored. This suggest that their emissions 
emerge from gas with the same physical conditions. The bottom 
plot shows the overlap zone. There is a large zone of excitation 
conditions where the column density of HNC is between 3 and 
10 times larger than that of HCN. The contour lines of -0.5 de- 
pict the zone of the physical conditions for which the column 
density of HNC is about 3 times larger than that of HCN. For 
both molecules and transition lines the optical depth ranges be- 
tween 0.001 and 30. The optically thick limit is basically defined 
by the Hmit of convergence of RADEX (Tmax ~ 100). 



According to the recent PDR and XDR models by Meijerink 
& Spaans (|2^5] | and Meijerink et al. (|2007] |. under PDR condi- 
tions a temperature = 80 K can be reached at a total column 
density of A^h ~ 3 x 10^' cm"^, if the Habing flux is Go ~ 10^ 
(a Habing flux Go = 1 corresponds to a far UV flux Ffjjy - 
1.6 X 10"^^ erg cm^^ s"'), or at A^h ~ 8 x 10^' cm^^ if Go ~ 10^ 



At those depths the density ^(Hi) would be slightly higher than 
10^ cm"^, though. On the other hand, in an XDR environment a 
temperature of 80 K and a density of 10^ cm"^ can be reached 
at a slightly larger column density of A^h ~ 10^^ cm"^, if the im- 
pinging radiation field is Ffuv ~ 1.6 erg cm"^ s"', or at a much 
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larger depth equivalent to A^h ~ 2 x 10^"* cm"^ if the radiation 
flux is Ffuv ~ 160 erg cm"^ s"'. 

However, at at a temperature of 80 K and density of 
10^ cm"^ there is no solution for HNC. This means that either 
the HNC emission arises from a gas with different physical con- 
ditions than HCN, or that both molecules trace a cooler (T^ < 80 
K) and lower density gas (n(H2) < 10^ cm"-'), if the emission of 
both molecules arise from the same gas. 

If the HCN and HNC emissions actually trace the same gas, 
the range of possible excitation conditions can be constrained 
to the zone where the conditions for HCN and HNC overlap. 
The bottom plot of Figure |6] shows the average overlap zone, 
where the contours correspond to the ratio between the total 
column densities of HCN and HNC. The overlap zone is de- 
fined by the optically thin limit (right edge of the excitation 
maps) of the HNC transition lines and the optically thick limit 
(left edge of the excitation maps) of the HCN lines. The mod- 
els show that A^(HCN) is larger than A^(HNC) in all the overlap 
zone. According to Figure 10 in Meijerink & Spaans (120051) . 
A^(HNC)/N(HCN) column density ratios lower than unity can 
be found mostly in PDRs, but also in XDR environments if the 
total column density A^(H) is lower than 10^"^ cm"^. If we still 
assume the same temperature (80 K) proposed by Tacconi et 
al. ( 1 19941 1, and we consider the contour line where A^(HCN) is 
about 2.5 orders of magnitude larger than A^(HNC), we find that 
the gas density should be n(H2) ~ 6 x lO'* cm"^, which agrees 
with the observed HNC 3-2/1-0 line ratio discussed in section 
4.2. However, this density would not be consistent with the re- 
sult found by Tacconi et al. (1994) nor with the observed HCN 
3-2/1-0 line ratio, which implies that the HNC emissions arise 
from a more diffuse gas than HCN. 



4.3.2. NGC 3079 

In contrast to NGC 1068, the excitation conditions modelled for 
the HCN and HNC line ratios observed in NGC 3079, overlap 
in most of the range explored (Figure |7]). This suggest that the 
emission from both molecules likely arise from the same gas. 
The spectral line shapes of HCN and HNC, showed in Figure [3] 
also hint that their distribution may be the same, although their 
line centers and widths are affected by noise. The bottom plot of 
Figure |7] shows the overlap zone of the excitation conditions. In 
most of these conditions A^(HNC) is between 0.5 and 1 order of 
magnitude larger than A^(HCN). 

According with the PDR and XDR models by Meijerink & 
Spaans (2005), a column density ratio A^(HNC)/N(HCN) larger 
than unity can be found only in XDR environments at a total 
column density A^(H) > lO^** cm"^, with radiation fields of 
about 1.6 erg cm"^ s"' (or equivalent Habing flux Go ~ lO-*) 
and total density n(H) ~ 10^ cm ^ A?(HNC)/N(HCN) column 
density ratios larger than unity at A^(H) > 10^"* cm"^ can also 



be found with stronger radiation fields (F^ ~ 160 erg cm s ) 
and higher densities (n(H) ~ 3 x 10^ cm"^^). On the other hand, 
A^(HNC)/N(HCN) can be larger than unity at N(H) < lO^^ cm^^ 
if the radiation fields are on the order of 1 .6 erg cm"^ s ' and the 
total density is about 3x10^ cm"-*. 

In our models we did not explore densities lower than 
10"* cm"^ since the HCN and HNC molecules are expected to 
trace higher densities. If we assume a kinetic temperature of 80 
K, as in the case of NGC 1068, and a density n(H2) ~ 10^ cm \ 
the column densities per line width would be about 10'^ and 
10'^ cm"^ km ' s for HCN and HNC, respectively. In the case 
of HNC both transition lines are equally optically thick at 80 K 
in the whole range of densities, whereas HCN is more optically 




1 ',! 




Fig. 8. Excitation conditions modelled for the 3-2/1-0 line ratio 
of HNC (top) and HCN (bottom) observed in NGC 2623 and 
NGC 1365, respectively. For NGC 2623 the optical depth of 
HNC 3-2 ranges between 0.003 and 10 in the J=l-0 line, and 
between 0.32 and 30 in the J=3-2 line. In the case of NGC 1365 
the optical depth in the HCN 3-2 line ranges between 0.1 and 
30, and between 0.01 and 10 in the J=l-0 line. 



thin. At a density of 10^ cm the optical depth of HNC is ~1 in 
both lines, whereas t is ~0.3 and ~0.1 in the 7=1-0 and J-3- 
2 lines of HCN, respectively. A density lower than 10^ cm"^ is 
more likely to be the right case in order to be consistent with the 
line ratios observed and discussed in section 4.2. 

4.3.3. NGC 2623 and NGC 1365 

Both NGC 2623 and NGC 1365 have similar excitation condi- 
tions for HNC and HCN, respectively. As in the case of HCN 
in NGC 1068, both molecules could be embedded in gas with 
densities larger than lO'' cm"^ at 80 K, which is consistent with 
line ratios larger than 0.3, as discussed in section 4.2. The col- 
umn densities per line width of these molecules would be about 
6 X lO'-' cm"^ km"' s for HNC and about lO'"* cm"^ km"' s for 
HCN, if the gas density is 10^ cm""*. 



4.4. The CN/HNC and CN/HCN line ratios 

The CN/HNC ratio increases towards the CND in both galax- 
ies, NGC 1068 and NGC 3079. The CN/HNC 1-0 ratio is lower 
than unity in all the galaxies, with the exception of NGC 1068. 
This can be interpreted as that HNC is more abundant than 
CN, assuming the same excitation conditions and that the emis- 
sion emerge from the same gas and distribution in the galaxy. 
However, since we do not really know the source size of CN (and 
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basically neither that of HNC) the beam dilution effects could be 
more (or less) severe than estimated here. 

Although the beam dilution uncertainty of CN also applies to 
the CN/HCN ratio, this case is more interesting because none of 
the galaxies where we do have data show ratios larger than one, 
neither for the global ratio (7= 1-0) nor for the nuclear part of the 
galaxies (higher transitions). All of these galaxies are considered 
active galaxies, so the presence of an AGN increases the chances 
of finding an X-ray dominated region (XDR) in their nuclear 
regions, as it seems to be the case of NGC 1068, according to 
Usero et al. (2004). 

In an XDR the CN molecule is expected to be more abundant 
than H CN (e.g., Lepp & Dalga rno 119961 Meijerink & Spaans 
12005] and Meijerink et al. l2007l) . and hence the CN/HCN inten- 
sity ratio could also be larger than unity, especially for the higher 
transition lines - CN 2-1 and HCN 3-2 - where the beam dilu- 
tion is less important than for the J-l-Q line. However, this is 
not seen in any of the sources presented here, nor in the sample 
of galaxies shown by Aalto et al. (2002), where results for AGN 
are presented along with starburst galaxies. 

According to Meijerink et al. (2007), in a PDR environ- 
ment the A^(CN)/A^(HCN) column density ratio ranges between 
2.0 and 0.5 for densities between lO** and lO*" cm"^, respec- 
tively, whereas in an XDR this ratio varies from over a 1000 (at 
n(H2) ~ 10^ cm-3) to 40 (at «(H2) ~ 10^ cm^^). If the CN inten- 
sity lines are proportional to the column density, and the estimate 
of the beam dilution is accurate enough, our results would favor 
a PDR scenario, rather than an XDR, with densities < lO*" cm"^^ 
in the central regions of all these galaxies. However, if the CN 
molecule is indeed more abundant than HCN, then the weaker 
intensity lines could be also due to a stronger optical depth ef- 
fect (in the escape probability sense) for CN. A rigorous model- 
ing and analysis of the CN molecule would be required in order 
to understand and predict the intensities of this molecule. 

On the other hand, the PDR scenario would also be consis- 
tent with the results found for the HNC and HCN molecules de- 
scribed above. HCN/HNC line intensity ratios larger than unity 
are usually found in PDR environments at total column densities 
A^(H) lower than 10^^ cm"^, while this ratio is larger for strong 
(> 10 erg s"' cm"^) radiation fields and low (~10"* cm""*) densi- 
ties in an XDR environment (Meijerink et al. l2007t . The lower 
densities required to observe this ratios in an XDR environment 
tend to dismiss this alternative since our simulations in section 
4.3 favor densities ~10^ cm"^ for the HNC and HCN molecules, 
as in the case of NGC 1068. In particular, there is evidence of 
recent starburst in the dense nuclear disk of NGC 1068 (Davies 
et al. 120061) ■ These results are consistent with a model consid- 
ering both, AGN and starburst components, required for mod- 
ellin g the U V to FIR atomic spectrum of NGC 1068 (Spinoglio 
et al. l2005T l. 

5. Conclusions 

We have used the SEST and JCMT telescopes to carry out a sur- 
vey of CN 2-1, HCN 3-2 and HNC 3-2 line emission in a sam- 
ple of 4 Seyfert galaxies, plus NGC 3079 which was observed 
with the IRAM 30m telescope. The conclusions we draw are as 
follows: 

1) We detected HNC 3-2 emission in 3 of the 5 galaxies, while 
we obtain an upper limit for one of them (NGC 7469). HCN 
3-2 was also detected in 3 galaxies (NGC 3079, NGC 1068 
and NGC 1365), while it was not detected in NGC 2623. CN 
2-1, along with the spingroups (J = 5/2 - 3/2, F = 7/2 - 



5/2) and {J = 3/2 - 1/2, F ^ 5/2 - 3/2) was also detected in 
NGC 3079, NGC 1068 and NGC 1365. 

2) The Hne shapes observed in NGC 1365 and NGC 3079 sug- 
gests that there is no circumnuclear disk in these galaxies. 

3) We find that in 3 of the galaxies the HNC 3-2/1-0 line ratios 
suggest that the HNC emissions emerge from gas of densi- 
ties n < 10^ cm"-^, where the chemistry is dominated by ion- 
neutral reactions. In NGC 2623 a model of large masses of 
hidden cold (10 K) gas and dust, as well as a chemistry dom- 
inated by ion-neutral reactions, are yet to be distinguished as 
the correct interpretation for the bright HNC observed in this 
galaxy. 

4) The 3-2/1-0 line ratios and the modelled excitation condi- 
tions imply that the HNC emission emerges from a more 
diffiise (n < 10^ cm"^) gas region than the HCN emission 
(n > lO'' cm"^) in NGC 1068, whereas they emerge from the 
same lower density (n < 10^ cm""*) gas in NGC 3079. 

5) The HCN/HNC and CN/HCN line ratios tentatively fa- 
vor a PDR scenario, rather than an XDR one, in the 3 
Seyfert galaxies where we have CN, HNC and HCN data. 
The A^(HNC)/A^(HCN) column density ratios obtained for 
NGC 3079 can be found only in XDR environments. 

In order to complete the sample, we plan to observe HCN 
3-2 and CN 2-1 in NGC 7469, CN 2-1 in NGC 2623 and HNC 
3-2 in NGC 1365. We plan to perform high resolution observa- 
tions to further study the distribution and source sizes of CN and 
HNC. 

Modeling of the collision data for the CN molecule would 
be useful to estimate the A^(CN)/A^(HCN) column density ratio, 
which would complement the A^(HNC)/A^(HCN) ratio in order 
to have a more sophisticated tool to estimate and distinguish the 
prevalent environment conditions of the high density gas in the 
nuclear region of Seyfert galaxies. 

The AGN contribution (through XDR effects) is typically of 
a small angular scale and can be seriously affected by beam di- 
lution at the transition lines studied in this work. On the other 
hand, the starburst contribution is of a larger angular scale than 
the AGN, and it effects can be contaminating our observations, 
and hence leading to the favored PDR scenario found with our 
models. Hence, our suggested interpretations could change if we 
zoom in on these sources. Therefore, high resolution maps of 
HNC and CN molecules are necessary to complement those of 
HCN, and to do a more accurate estimate of molecular abun- 
dances and line intensity ratios, which take source size into ac- 
count. Observations of the higher transition lines (e.g. 7=4-3) 
can also aid to disentangle the effects of the AGN and the star- 
burst ring, due to the smaller beam size obtained at higher fre- 
quencies. 
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